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Vlie  robustness  of  stability  conditions  for  linear 
time-invariant  feedback  systems  is  examined 
assuming  three  different  types  of  representations; 
state«>space  representation,  coprime  matrix 
fraction  representation,  and  transfer  Juj^ctj-cyi^ 
representation.  otres^  the  importance^of 

certain  details  of  the  representation  used  and, 
even  more,  the  importance  of  making  sure  that 
the  allowed  perturbations  be  relevant  to  the 
physical  situation  under  study. 

*v 

1.  Introduction 

Engineers  design  for  production:  therefore  it 
is  required  that  their  nominal  design  as  well  as 
a very  high  portion  of  the  systems  produced  — 
which  suffer  from  element  deviations  and  manu- 
facturing tolerances  — meet  the  specifications. 
Furthermore,  the  systems  produced  must  meet  the 
specif ications  not  only  as  they  leave  the 
production  line  but  also  in  the  field  — where 
they  suffer  from  temperature  effects,  aging, 
weathering  etc. — . Hence  the  interest  in 
sensitivity  and  robustness.  It  is  for  these 
reasons  that  this  subject  has  an  extensive 
literature  [e.g.  1,2]. 


relevance  of  some  parameter  perturbations. 

In  order  to  avoid  repetitions,  we  start  by  defining 
some  terms  and  some  notations. 

II.  Preliminary  Definitions 

IR  , C,  IR  (s)  and  IRfs]  denote,  respectively,  the 
fields  of  real  numbers,  of  complex  numbers,  of 
rational  functions  with  real  coefficients  and 
the  commutative  ring  of  polynomials  with  real 
coefficients.  The  superscripts  "n"  and  "nxm" 

(as  in  IR  , IR  (s)  ) denote  the  corresponding 

ordered  n-tuples  and  nxm  arrays. 

« {s  ^ C I Re  s ^ 0}  denotes  the  closed  right- 
half-plane.  = (s  ^ (E|Re  s < 0}  denotes 

the  open  left-half-plane.  Given  any  scalar 
rational  function,  we  assume  once  and  for  all  that 
it  is  written  as  n(s)/d(s)  where  the  polynomials 
n and  d are  coprime  and  d is  nonic, 

A continuous-time,  linear,  time-invariant,  lumped, 
multi-input  multi-output  system  is  said  to  be 
exponentially  stable  (abbr.  exp.  stable)  iff  its 
transfer  function  C(s)  G IR  is  proper 

(i.e.  bounded  at  infinity)  and  G(s)  has  no 
I^-poles.  For  example,  for  the  system  shown  on 
Fig.  1,  this  means  that  G;  (uj^,U2)' — ►(e^,e2)  has 
these  properties. 


In  this  paper,  we  consider  the  robustness  of  the 
stability  conditions  for  a continuous-time, 
linear,  time-invariant,  lumped,  multi-input 
multi-output  feedback  system.  (See  Fig.  1)  If 
the  feedback  system  is  made  of  an  interconnection 
of  stable  subsystems,  it  seems  intuitively  clear 
that  under  some  reasonable  conditions  and  under 
reasonable  allowed  perturbations  the  stability 
conditions  are  robust.  But  what  if  the  subsystems 
are  unstable?  Might  it  not  happen  that  due  to 
perturbation  some  kind  of  pole-zero  cancellation 
is  destroyed?  The  purpose  of  this  paper  is  to 
examine  the  conditions  vjhich,  under  several 
representations,  guarantee  robustness  of  the 
stability  conditions.  We  will  find  that  the 
nature  of  the  representation  and  the  nature  of 
the  allowed  perturbations  play  a crucial  role. 

This  will  also  lead  us  to  make  some  remarks  on  the 

+ 
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III.  System  Description 


We  consider  the  input-output  stability  problem  of 
the  continuous-time,  linear,  time-invariant,  lumped, 
multi-input  multi-output  feedback  system 
S:  (u.  ,U2)» — >(e.  ,e~)  described  in  frequency  domain 
by  (see  Fig.  1) 


^ V2 


- Vi 


where  C^,  C2  ^ lR(s)"’‘",  the  u^'s  are  the  inputs 
and  the  ej^'s  are  the  "errors".  The  transfer  function 
of  S Is  G:  (uj^.u,)! — ►(e,.e,),+  Throughout  this 
this  paper,  we  mate  the  following  assumption: 

Assumption:  The  transfer  functions  >^2  proper 
(i.e.  bounded  at  infinity)  and 


We  only  need  to  consider  (u^,U2)^ 


.(Cj^.e^)  because 


the  map  (Uj .u^).— *(y ,y2)  is  exp.  stable  if  and 
only  if  the  map  (u  ,u  )t — »(e, ,e  ) is  exp.  stable 
(lOJ.  ^ ^ ^ ^ 


Appcofd  for  public  rel®a»«: 
Dlitilbutloi)  Unlijnitwi 


det[I-H:2(")  0.  (2) 

This  assumption  Implies  that  the  transfer  function 
G of  S exists  and  Is  proper.  Note  that  a transfer 
function  must  be  proper  In  order  to  have  a 
state-space  representation. 

We  Investigate  the  following  question:  given  that 
Gj,G2  are  described  In  a specified  way  and  given 
some  set  of  allowed  perturbations,  are  the 
Input-output  stability  conditions  for  the  system 
S robust  under  the  allowed  perturbations  in 
Gj^,G2?  By  this  we  mean:  if  the  feedback  system 
S is  exp.  stable  for  the  nominal  values  of 
G^,G2,  does  It  Imply  that  it  will  remain  exp. 
stable  for  all  sufficiently  small  perturbations 
in  Gj^,G2,  selected  from  the  allowed  set?  It  will 
turn  out  that  the  answer  depends  very  much  on 
the  representation  of  Gj^,  G2  and  the  set  of 
allowed  perturbations. 

IV.  A Preliminary  Lemma 

All  the  argumentation  below  hinges  on  a lemma  that 
essentially  says  "small  perturbations  in  the 
coefficients  of  an  algebraic  equation  cause 
small  perturbations  in  its  zeros."  More  precisely, 
we  state  a well-known  lemma. 

Lemma  [3,  Thm  9.17.4].  Let  D(z^;e)  denotes  the 
open  disc  in  C centered  on  and  with  radius  e. 
Consider  the  polynomial  p defined  by 

p(s)  “ an®"  ^ + ...  + a ,s  + a (3) 

^0  1 n-i  n 

where  a^  6 IR  , Vi,  and  without  loss  of  generality, 
Oq  > 0.  Let  the  polynomial  p have  q pairwise 

distinct  zeros:  zj,Z2, ,z^  with  respective 

multiplicities  m^^,  (hence  = n) . Then  for 
all  e > 0,  there  is  an  n(e)  > 0 such  that  for 
all  Saj.  satisfying 

jSaj^l  < n i “ 0,1 n (4) 

the  perturbed  polynomial 

(an+fSa^ls^  + ...  + (a  ,+6a  ,)s  + (a  +6a  ) 

'00  n-1  n-1  n n 


still  has  m^  zeros  in  D(Zj^;e)  for  i 


1,2,... ,q. 


Comments:  (1)  If  all  zeros  of  p were  simple, 
this  Lemma  would  be  a direct  consequence  of 
the  Implicit  function  theorem;  the  point  is  that 
the  continuous  dependence  of  the  zeros  is  still 
valid  in  the  case  of  multiple  zeros.  (il)  It  is 
crucial  to  observe  that  the  degree  of  p was  not 
affected  by  the  perturbations:  Indeed  suppose 
that  Instead  of  (5)  we  had 


6a  ^ 


+ (a„+6a^) 


as  the  perturbed  polynomial;  then,  for  n > 0 
sufficiently  small,  if  |6aj|  < n for 
I ■ - 1 , 0, 1 , • • . ,n,  p would  nave  n-fl  zeros,  n of 
them  in  the  discs  D(z£ic)  and  one  approximately 
equal  to  -(aQ+6aQ)/6a_j.  (This  approximate  zero 
is  the  leading  term  of  a sequence  of  successive 
approximation  which  converges  for  n small). 


This  additional  zero  is  very  large:  its  sign  is 
positive  or  negative  according  to  6a^j^  < 0 or 

> 0,  respectively.  Similarly,  if  there  had 
been  several  additional  terms  of  degrees  larger 
than  n,  there  would  have  been  several  such  zeros 
with  very  large  magnitude  and  whose  location  in 
the  s-plane  depends  on  the  magnitudes  and  signs 
of  the  6a_j^*s.  Note  that  when  the  perturbed 
problem  is  of  the  form  (6),  we  are  essentially 
dealing  with  a singular  perturbation  problem; 
see  e.g.  [4,5 j . 

V.  Robustness  Results 

Case  I : The  transfer  functiorisG|  and  G-^  are 
specified  by  minimal  state-space 
representations 

It  is  well  known  that  if  [A,B,C,D]  is  a minimal 
state-space  representation  of  the  proper  transfer 
function  G,  then  G is  exp.  stable  if  and  only 
if  det(sI-A)  has  all  its  zeros  in  l_.  For 
i = 1,2,  let  [A^,Bj^,Cj,Dj^J  be  a minimal  state- 
space  representation  of  the  proper  transfer 
function  with  the  state  x^^  ^ IR^^.  Let 
[A,B,C,D]  be  the  state-space  representation  of  the 


feedback  system  S:  (u^,U2)»- 


with  the 


state  (Xj^jX^).  We  know  that  [6]  (i)  [A,B,-,D]  is 
a minimal  state-space  representation  of  S,(ii) 


B2  ( I+D^D^)  A^-B^  ( I+Dj^D^) 


and  (Hi) 

det(sI-A)  = det(sI-A^) 'detCsI-A^) ’detCl+G^G^) (s) 


Now  the  feedback  system  S is  exp.  stable  if  and 
only  if  det(sI-A)  has  all  its  zeros  in  C_; 
furthermore  for  any  perturbation  [6Ai,6Bj^,6Ci,6DiI 
in  the  constant  matrices  lAi,Bi,C£,Di] , i * 1,2, 
the  degree  of  det(sI-A)  remains  equal  to 
In  view  of  (7),  we  obtain  from  the  lemma  above 
a well-known  result: 

Robustness  Result  1.1. 

If  the  given  feedback  system  S is  exp.  stable  at 
the  nominal  data  point,  then  for  any  sufficiently 
small  perturbation  [6Ai,6Bi,6Cj^,6Di] , i * 1,2, 
the  resulting  perturbed  system  is  still  exp.  stable. 

□ 

Remark  I . 1 : In  many  applications,  neither  measure- 
ments nor  system-component  element  values  directly 
specify  the  [Aj^,Bi,C^,D^] ; therefore,  engineers 
should  ask  whether  the  perturbations  [6A^,6Bj^, 60^,60^1 
of  this  analysis  cover  all  the  possible  perturbations 
expected  in  the  contemplated  physical  environment. 

Remark  1.2:  Note  that  for  sufficiently  small 
allowed  perturbations,  [ Ai+6Aj , B j^+6Bj ,C^+6C , 

Di+6Di]  will  also  be  a minimal  state-space 


♦ 
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representation.  However  Robustness  Result  I.l 
still  holds  even  if  B£+6B^,Cj^+6C^,Dj^+6Dj^ ) 

Is  not  minimal. 

Case  II:  The  transfer  functlonsG^  and  G„  are 

specified  by  coprime  matrix  fractions 
representations 

Let  be  specified  by  their  coprime 

factorizations: 


G,  - N,  D, 

1 ir  Ir 


where  Nj^^,  Dj^^,  ^ K [s]  and  the 

pairs  are  rlght-coprirae  and 

lefc-coprime  respectively  [7,8,9].  Let  the 
allowed  parameter  perturbations  be  perturbations 
in  the  coefficients  of  each  scalar  polynomial 
entry  in  the  four  polynomial  matrices 

^2Jl’  without  increasing  the  degree  of  any 
scalar  polynomial. We  know  that  [10,11]  the 


feedback  system  S:  (u^jU^)^' 


is  exp. 


stable  if  and  only  if  its  characteristic 
polynomial 

a:  = det[D2^D^^  + (9) 

has  all  Its  zeros  In  5_.  We  can  also  write 
A = det  D^^-det  D^_^-det[I  + d' ] 

= det  D^^'det  Dj^^'det[I  + (10) 

By  assumption  (2),  G^,  are  proper  and,  for 
large  s,  det(I  + G2Gj^](s)  - det[I  + G2(®)G^(®)] 

+ 0(l/s).  Hence  degree  of  A = degree  of 

det  + degree  of  det  0^^.  In  view  of  (10)  it 

follows  from  the  lemma: 


D(s) , Then  degree  of  det  D £ min  { r , , ^ c , 

i»l  ^ 1“1  ^ 

n 

If  the  equality  holds  with  ^ r.t  then  D Is  said 

i-l  ^ 

to  be  row-proper  and  If  the  equality  holds  with 
n 

y!  c . , then  D is  said  to  be  column-proper  [7,8]. 
i*l 

Thus  if  D is  either  row-proper,  or  column-proper, 
then  it  is  easy  to  see  that  the  degree  of  det  D 
will  never  be  increased  by  any  allowed  perturbation 
6D.  It  is  well-known  that  [7,8]  for  any  nxn 
polynomial  matrix  D(s)  with  det  D(s)  t 0,  there 
exist  unimodular  matrices  U^(s),  U^(s)  such  that 
D(s)  U^(s)  is  column-proper  and  U^(s)  D(s)  is  row- 
proper.  Therefore  given  a rational  matrix,  there 
exist  a left-coprime  factorization  (N^,Dj^)  such 
that  is  row-proper  and  a right-coprlme 
factorization  (N^,D^)  such  that  is  column-proper. 
Summarizing  the  facts  above,  we  have 

Robustness  Result  II. 2. 

Suppose  that  row-proper  and  is  column- 

proper;  if  the  given  feedback  system  S is  exp. 
stable  at  the  nominal  data  point,  then  for  any 
sufficiently  small  allowed  perturbation 
6D.,^,6N2j^)  the  resulting  perturbed  system  is  still 
exp.  stable.  ^ 

Remark  11,2:  The  requirement  that  D...  and  D.  be 

T 21  Ir 

resp.  row-proper  and  column-proper  is  very 
important,  for  otherwise  the  degree  of  the 
characteristic  polynomial  may  increase  as  a result 
of  an  arbitrarily  small  allowed  perturbation. 

This  is  shown  in  the  following  example. 


Robustness  Result  II. 1. 

If  the  given  feedback  system  S is  exp.  stable  at 
the  nominal  data  point,  then  for  any  sufficiently 
small  allowed  perturbation  ('^^^r’^^lr 
which  does  not  increase  the  degree  of  det  D2jj^  and 
det  the  resulting  perturbed  system  Is  still 

exp.  stable.  ^ 

Remark  II.  1:  Note  that  for  sufficiently  small 
allowed  perturbations,  (N]^^+6N^^,Dj^j.+6D.  ) will 
be  right-coprime  and 

left-coprime.  However  Robustness  Result  II. 1 
still  holds  even  if  they  are  not  coprime. 


Example : Let  D« 


s+2  s+1 

s+1  s+3^ 


Hence  A(s)  = 1^2)1  ^lr^^2Jl^lr  ^ ” 5s+ll.  Now 

is  not  column-proper:  consider  a small 
perturbation  say,  a in  tne  coefficient  of  s in  one 


of  the  diagonal  elements  of  D, 


For  example. 


when  the  (1,1)  element  becomes  (l+a)s  + 2 we  obtain 

A(s)  + 6A(s)  = as^  + (5+Aa)s  + 11 

which  has  a zero  with  positive  real  part  whenever 
a < 0. 


Let  r^  (resp.  c^)  be  the  highest  power  of  s in  the 
^th  row  (resp.  column)  of  an  nxn  polynomial  matrix 

4>4- 

viewing  the  matrix  entries  as  polynomials,  note 
that  the  zero  polynomial  has  degree  and  a non- 
zero constant  polynomial  has  degree  0.  Thus  the 
allowed  perturbations  will  not  cause  a zero  in 
any  of  the  four  polynomial  matrices  to  become 


Case  HI:  The  transfer  functions  G^^  and  G^  are 
specified  bv  the  2n^  scalar  rational 
functions  which  are  the  elements  of 
and  G^. 

Let  the  allowed  parameter  perturbations  be  the 
perturbations  in  the  coefficients  of  the  numerator- 
and  denominator-polynomials  of  each  scalar 
rational  function  entry  in  the  matrices  Gj^,  G^ 
subject  to  the  condition  that  they 


White  Secthon 
Buff  Section 


DISTRIB!lTION/AYAyBlllTT  CODES 

iPist.  •■  All,  and/or  SPECIAI 


□ □ 


do  not  Increase  the  decree  of  any  scalar 
polynomial. 

I I I . 1 ; Slngle~lnput  single-output  subsystems 

Consider  the  case  where  » n^/d^^,  i » 1,2, 
where  n^,  d^  are  coprime  scalar  polynomials,  i.e., 
they  are  coprime  elements  of  IR[s].  Clearly  the 
characteristic  polynomial  of  the  feedback  system 
S is  (dj^d^+nj^n^) . Since  d^,  d^  are  scalar 
polynomials,  they  are  both  row-proper  and  column- 
proper.  Thus  by  Robustness  Result  11.2,  if  the 
given  feedback  system  is  exp.  stable  at  the  nominal 
data  point,  for  any  sufficiently  small  allowed 
perturbation  (for  Case  III),  the  resulting 
perturbed  system  is  still  exp.  stable. 

III. 2:  Multi- input  multi-output  subsystems 

Recall  the  notations  introduced  in  Case  I and 

det(sI-A)  = det(sI-Aj^)‘det(sI-A2)*det(I+G2G^)  (s) 

(8) 


where  the  sums  are  taken  over  all  the  (E^-poles 
X of  Gj^  and  G^,  or  equivalently,  over  all  the 
poles  X in  C^  of  G^  and  G^.  With  this  in  mind  we 
see  that  some  sufficiently  small  allowed  perturbation 
may  change  the  required  number  of  encirclements  in 
only  one  way:  namely  to  have 

P„+  T.  E6(Gi+4G  :A) 
i-1,2  * ^ ^ 

where  the  sums  are  taken  over  all  the  poles  X in 
C^  of  perturbed  transfer  functions  + 6G^ 
and  G2  ■*"  *^^2* 


Fur  example,  consider 


G^(s) 


j_p  ‘ 

s-ll_3  1.5 


, G^(s)  + 6Gj^(s)  • ^ 


•2+a  1 


then  A(Gj^;l) 


rank 


3 1.5. 
1,  [15, p. 115],  but. 


for  any  a ^ 0,  A(Gj^+6Gj^:  1)  - 2. 

Therefore  we  must  formulate  our  result  as  follows: 


Let 

p = number  of  C -zeros  of  det(sl-A, ) 

0+  + 1 

•det(sI-A2),  counting  multiplicities. 

(11) 

Then  recall  the  Graphical  Stability  Conditions: 
[12,13]  the  feedback  system  S:  (Uj^,U2)»— *(ej^,e2) 
is  exp.  stable  if  and  only  if  the  Nyquist 
diagram  of  s^— *det[I  + G2(s)G^(s)J  — for  the 
contour  C which  is  duly  indented  to  the  left  at 
all  jw-axis  poles  of  det[l  + G2(s)Gj^(s)]  — does 
not  go  through  the  origin  and  does  encircle  the 
origin  p^^  times  in  tlie  counterclockwise  sense. 

Now  suppose  that  for  the  nominal  parameter  values 
in  and  G2  the  Nyquist  diagram  satisfies  Che 
stability  conditions  above.  Consider  the  effect 
on  the  Nyquist  diagram  of  small  allowed  parameter 
perturbations  in  G^^  and  G2.  Now,  (a)  for  each 
s ^ (C,  except  at  poles,  det[I  + G2(s)Gj^(s)]  is  a 
continuous  function  of  all  the  numerator-  and 
denominator-coefficients  of  G^^  and  G2;  (b)  the 
Nyquist  diagram  of  the  contour  C is  a compact 
curve  in  t,  therefore  for  any  sufficiently  small 
allowed  parameter  perturbation,  the  Nyquist 
diagram  will  still  avoid  the  origin  and  encircle 
it  p^_j^  times;  (c)  for  sufficiently  small  allowed 
perturbations  and  for  a fixed  contour  C (indented 
Co  the  left),  any  jw-axis  pole  of  G^  and/or  G2 
chat  is  perturbed  will  still  remain  in  the 
open  set  in  C enclosed  by  the  indented  contour  C. 

Recall  that  Che  order  of  zero,  X,  of  dec(sI-A^) 
is  equal  to  Che  McMillan  degree  of  G^  at  its 
pole  X — which  we  denote  by  A(G. ;X)  [14,15]. 

Thus  ^ 

p„.  ” E E 

° 1-1,2  A ^ 


Robustness  Result  III.l. 

If  the  given  feedback  system  S is  exp.  stable  at 
the  nominal  data  point,  then  for  any  sufficiently 
small  allowed  perturbations  6G^,  6G2  which  satisfy 

E E^(g,;a)=  E E ‘''(g^+6g.:A) 

1=1,2  X 1=1,2  X 

where  the  sums  are  taken  over  all  the  poles  X in 
of  the  corresponding  transfer  functions,  the 
resulting  perturbed  system  is  still  exp,  stable. 

n 

In  parcicularf  robustness  of  stability  conditions 
for  the  feedback  system  S follows  if  and  G2  are 
both  exp.  stable. 

In  the  following  discussion,  we  restrict  ourselves 
sftnple  C^-poles  of  G^^  and  G2  because  there 
always  exists  some  allowed  perturbation  which 
splits  a multiple  pole  of  G^  (i=l  and/or  2)  in 
many  ways.  To  see  this  suppose  Chat  the  (1,1) 
element  of  has  a third-order  pole  at  s=  -1,  thus 
its  denominator,  d^^,  has  the  form 

dii(s)  = (s+1)^  Pll^®^  ^11^^^  ^ ^ 

This  triple  pole  can  be  split  into  three  simple 
poles,  for  example 

d^^(s)  + 6dj^j^(s)  - [(s+1)^  + 

or  into  a double  pole  and  a simple  pole 

dii(s)  + 5d^j^(s)  = (s+l)^[(s+l)  + ejpj^jCs), 

etc.  More  generally,  d^^j^  has  multiple  zero(s)  if 
and  only  if  its  discriminant  = 0 [16). 

Since  is  a polynominal  in  the  coefficients  of 
it  is  clear  that  multiple  zeros  are  not 
generic  [2]. 


Consider  now  the  effect  of  an  arbitrary  small 
allowed  perturbation  on  Che  McMillan  degree  of 
^1  ^k*  ^ simple  pole  with 

Usually  only  some  of  the  n^  elements 
of  G|  have  a pole  at  and  for  any  sufficiently 
small  allowed  perturbation,  if  the  (l,k)  element 
of  has  iw  pole  at  A^,  then  the  (l,k)  element 
of  Gj^+6G^  will  still  have  no  pole  in  a small 
neighborhood  of  A^^.  Consequently,  only  the 
nonzero  elements  in  the  residue  matrix  of 
Gi  at  are  affected  by  the  allowed  perturbations* 
Thus  R,  is  a structured  matrix  in  the  sense  of 
Shields  and  Pearson  [17,18]  i.e.,  it  has  a fixed 
pattern  of  zero  elements.  Let  be  the  number 
of  coefficients  which  specify  G^.  The  generic 
rank  [17,18]  of  the  structured  matrix  R^^  is 
def ined  to  be  the  maximal  rank  that  achieves 
as  a function  of  these  parameters.  does 

not  achieve  its  generic  rank  only  for  parameter 
values  in  some  proper,  closed,  nowhere-dense 
variety  V C IR^I.  Consequently,  if  the  rank  of 
the  nominal  R^^  is  less  than  its  generic  rank, 
then,  for  some  arbitrarily  small  allowed 
perturbation,  the  rank  of  R^^  will  jump  to  its 
generic  rank.  Thus  we  have  the 

Robustness  Result  III. 2. 


Suppose  that  all  C^-poles  of  G^^  and  of  are 
simple  and  that  for  each  of  these  poles  the 
nominal  residue  matrix  has  a rank  equal  to  its 
generic  rank;  if  the  given  feedback  system  is 
exp.  stable  at  the  nominal  data  point,  then  for 
any  sufficiently  small  allowed  perturbation,  the 
resulting  perturbed  system  is  still  exp.  stable. 

Furthermore,  if  the  nominal  residue  matrix  of  some 
C^-pole  of  either  G^^  or  G2,  has  a rank  less  than 
its  generic  rank,  and  if  the  nominal  system  S is 
exp.  stable,  then  for  some  arbitrarily  small 
allowed  perturbation,  the  resulting  perturbed 
system  is  unstable.  ^ 

Remark  IIl.l:  The  allowed  perturbations  considered 
in  Case  111  appear  quite  reasonable.  However,  we 
should  be  on  guard  that  they  might  include 
perturbations  that  have  no  physical  meaning  for  the 
case  at  hand.  For  example,  this  could  occur  if 
G. , instead  of  being  specified  by  a collection 
of  rational  functions  — the  entries  of  the 
matrix  — were  specified  by  a block  diagram 
delineating  the  interconnections  between  the 
subsystems  constituting  G^.  In  that  case  the 
appropriate  perturbations  to  consider  are  not 
arbitrary  perturbations  in  all  the  coefficients  in 
the  n^  rational  function  specifying  G.  but  rather 
perturbations  in  the  parameters  specifying  the 
subsystems  constituting  Gj^.  This  is  illustrated 
by  the  following  example. 

Example.  Suppose  that  G^  consists  of  a collection 
of  subsystems  in  series  and  in  parallel  (no 
feedback!)  and  that  only  one  subsystem  has  an 
unstable  pole,  say  at  p,  with  Re  p ^ 0.  The 
contribution  of  that  subsystem  to  Gj^  is  exhibited 


on  Fig.  2:  the  l.th  scalar  input  of  G can  only 
affect  the  scalar  input  v of  the  unstable  sub- 
system through  the  gain  similarly  the  unstable 
subsystem  output  z is  also  scalar  and  is  fanned  out 
to  the  ^th  output  of  Gj^  through  gain  y^.  Clearly 
the  residue  of  G^^  at  p is  the  dyad  (i.e.  a rank- 

one  matrix!)  R “ where  y - (Y,,yo,---,y  )^ 

X n 

B'  = (B^.B2,  • • • ,6jj)  • Now  any  small  perturbation 
in  the  physical  parameters  (y,,B,),  i • 1,2, ...,n 
and  p will  not  change  the  rank  of  R because  the 
dyadic  structure  of  R is  dictated  by  the  nature  of 
the  interconnection  of  the  subsystems  constituting 
Gj^,  and  not  by  the  numerical  value  of  the  parameters. 


Another  way  of  viewing  this  fact  is  to  say:  the 
only  meaningful  physical  parameters  which  determine 

R are  the  2n  scalars  y,,  y.,...,y  , B,,B.,...,B  . 

12  n 1 2 n 

To  take  the  abstract  mathematical  point  of  view 
that  R is  an  array  of  n^  real  numbers,  and  there- 
fore n^  independent  perturbations  of  its  parameters 
are  appropriate  is  mathematical  fiction,  and  not 
an  analysis  of  the  physical  system  under  consider- 
ation. 


The  lesson  of  this  remark  is  that  before  consider- 
ing robustness,  one  should  go  back  to  the  physical 
model  of  Che  system  under  consideration  and  trace 
out  the  effect  of  perturbations  of  the  physical 
parameters  of  the  model  on  the  coefficients  of  the 
mathematical  representation.  Only  in  this  way, 
the  engineer  will  assure  himself  that  the  allowed 
perturbations  he  worries  about  pertain  to  physical 
reality  and  not  to  mathematical  fiction. 


Remark  III. 2:  There  is  another  case  where 
perturbations  are  more  restricted  than  those 
considered  in  Case  III:  in  many  models  Newton's 
law  dictates  a second  order  pole  at  the  origin  — 
say  in  the  transfer  function  from  the  external 
forces  and  the  center  of  mass  — . Clearly  this 
second  order  pole  is  not  subject  to  perturbations 
due  Co  measurement  or  manufacturing  errors! 


Example:  Consider  the  planar  mechanical  system 
described  by 


myi  - ky^  = Uj^ 


my 


"2- 


where  (Uj^,U2)  is  the  applied  force,  (yj^,y2)  is  the 
position  of  the  particle  of  mass  m and  the 
physical  parameters  are  m > 0 and  k > 0.  Here 


“ -1  -2 
m s 

-2. 

m ks 

-2 

as 

Bs‘^ 

0 

-1  -2 
m s 

0 

as-i 

where  we  put  a - m"^,  6 • m“^k.  Clearly,  in  this 
example  the  second  order  pole  of  G and  the  fourth 
order  pole  of  G are  not  subject  to  perturbations 
when  the  physical  parameters  a,  6 are  perturbed. 

In  the  present  case,  the  McMillan  degree  at  the 
unstable  pole  (s“0)  is  insensitive  to  small 
perturbations  in  the  physical  parameters.  To 

« k-4 

check  this  write  G(s)  ■ V R.  s ; then  A(G;0)  is 
k-0 

given  by  the  rank  of  the  Hankel  matrix  [151 
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H = 


o 


■'3  ■'2  1 o 

R3  R^  0 

R,  R 0 0 

1 o 

R 0 0 0 


where  R 


12  S]  ■ ■ C 3 


and  the  others 


o _ 

are  zero.  It  is  easily  checked  that  for  all 
allowed  physical  values  (namely,  o > 0,  6 > 0) 
A(G;0)  = 4 and  also  that  the  generic  rank  of 
is  4,  (17,  p.  211].  Thus  if  this  transfer  funct 
G were  the  forward  gain  of  a closed-loop  system 
with  an  exponentially  stable  feedback  gain 
K(s)  - diag(k  (s),k2(s))  then  p^,  the  number  of 


required  counterclockwise  encirclements,  would  be 
4 and  it  would  be  robust  under  any  perturbation  of 
the  physical  parameters  o and  6.  j, 


Relation  between  the  effect  of  parameter 
perturbations  in  the  three  cases  above. 

Consider  a given  transfer  function  G(s),  one  of 
its  minimal  state-space  realization  [A,B,C,D]  and 
one  of  its  right-coprlme  factorizations  (N^,D^). 
Offhand  we  might  think  that  if  some  allowed 
parameter  perturbation  in  Case  I leads  to  the 
transfer  function  G+6G,  then  there  exists  some 
allowed  parameter  perturbation  in  Case  II  which 
leads  to  the  same  transfer  function  G+6G,  and 
so  on.  In  other  words,  we  might  think  that  the 
parameter  perturbations  in  the  three  cases  are, 
in  some  sense,  equivalent.  This  is  not  the  case: 
Indeed  Remark  II. 1 shows  that  for  some  coprime 

factorization  (N  ,D  ) there  is  some  aribtrarlly 
r r 

samll  allowed  parameter  perturbation  (of  Case  11) 
which  will  increase  the  degree  of  the  characteristic 
polynomial,  this,  however,  is  impossible  under 
any  allowed  parameter  perturbation  in  Case  I. 

Similarly  the  example  below  shows  that  for  some 
G(s)  there  is  some  arbitrarily  small  allowed 
parameter  perturbation  of  Case  I which  leads  to 
some  transfer  function  G+5G  unattainable  by  any 
allowed  parameter  perturbation  of  Case  III. 


Example:  Consider  a transfer  function  G(s)  with 
the  minimal  state-space  representation  [A,B,C,D) 
where 


A = 


Thus, 


1 

0 

-1 

0 

_-l 

G(s)  - 


0 0 0 0 

110-1 
0 0 0 0 

0 0 0 0 

0 0 0 0_ 

(sI-A)-l 


^5x5’  “ 


‘'5x5’ 


The  (1,2)  position  minor  of (sI-A)  is  given  by 


0 

0 

s 

0 


for  all  s. 


Therefore  the  numerator  polynomial  of  the  (2,1) 
element  of  G(s)  has  a degree  of  -“.  Let  SA  be  a 
5x5  matrix  which  has  only  one  nonzero  element  a 


at  (3,1)  position.  Clearly  [6A, 0,0,0]  is  an 
allowed  parameter  perturbation  in  Case  1.  The 
(1,2)  position  minor  of  (sl-A-6A)  is  given  by 


det 


"0  -10 
1-a  s 0 
0 Os 

Ll  0 0 

det[sI-(A+6A) ] 


1 

0 

0 

s 


Also 


Therefore  the  numerator  polynomial  of  (2,1)  element 
of  &f$G  has  degree  zero,  (Note  the  cancellation!) 
Due  to  the  increase  in  degree  of  the  numerator 
polynomial  of  (2,1)  element,  C+6C  is  unattainable 
from  G by  any  allowed  parameter  perturbation  of 
Case  III. 


VI.  Conclusion 


In  this  paper  we  studied  the  robustness  of  the 
exponential  stability  of  continuous-time,  linear, 
time-invariant,  lumped,  multi-input  multi-output 
feedback  systems.  We  presented  robustness  results 
for  three  types  of  system  representations  with 
corresponding  secs  of  allowed  parameter 
perturbations . 

The  robustness  results  above  were  obtained  assuming 
that  the  two  subsystems  G^  and  G.  have  the  same 
number  of  inputs  and  outputs,  all  the  results  can 
be  extended,  after  simple  modifications,  to  the 
case  where  the  number  of  inputs  and  outputs  of  each 
subsystems  are  different. 

Also  since  all  Che  arguments  used  are  purely 
algebraic  and  are  based  on  simple  properties  of 
rational  functions  and  polynomials,  all  the 
results  above  apply  equally  well  to  Che  discrete- 
time case  except  that  in  Robustness  Result  III.l, 

C should  be  interpreted  as  the  unit  circle  with 
inward  indentation  Co  avoid  the  poles  of  G^  and  G^, 
lying  on  the  unit  circle,  and  C^  should  be 
interpreted  as  the  "outside"  of  C (more  precisely 
Che  unbounded  connected  component  of  (t-C) . 
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